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the rates for -f*p — > 7X and for 7*7* — > 77, where X denotes the proton dissociation. We focus on 
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perturbation theory. However, our calculations do allow us to study the — t — > behaviour of the 
7*7* — > 77 process in the region where the incoming photons are sufficiently virtual. 
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1 Introduction 

One of the major challenges in the domain of strong interaction dynamics is to gain an understanding 
of diffr active phenomena within the framework of QCD. Typically, diffractive processes are driven by 
dynamics at low momentum transfers and hence cannot be studied using perturbation theory. 

Large-t diffraction has been shown to be one area where perturbation theory can be used |l]]. 
Double dissociation in hadron-hadron and photon-hadron collisions with a large momentum transfer 
between the diffracting systems has been, and continues to be, studied both theoretically [Q, |3| and 
experimentally j|]. The large momentum transfer means that jets are produced in the diffracting 
systems, the identification of which leads to a sample of 'gaps between jets' events. A reliable calcu- 
lation of the expected rate for double dissociation is made difficult because there is the possibility of 
additional strong interactions between the diffracting systems which could result in the loss of gaps. 
A process which avoids the 'gap survival' problem is diffractive vector meson production in photon- 
hadron collisions: 7£> — ► VX j|, ||, |?J ||. The meson V is produced with a high pt relative to the 
incoming photon and the largeness of this momentum transfer typically leads to the dissociation of 
the hadron into system X. A principal challenge in computing the rate for meson production lies in 
determining how the meson is produced. To avoid both the meson wavefunction problem and the gap 
survival problem, one can study the diffractive production of photons ||, ^, [|. At HERA, one can 
look for *yp — ► ^X where the photon is produced with a high pt and the proton dissociates. At LEP2 
or a future linear collider one can look for 77 — ► 7 A 7 " or 77 — ► 77 H, |l0| . In this paper we compute 
the rates for diffractive photon production in photon-proton and photon-photon collisions at large 
momentum transfer. Our calculations include the possibility of off-shell incoming photons. 



2 The 7*7 impact factor 

The scattering amplitude in the high energy (Regge) limit may be factorised in the following way 

A(s,t) = ^^k^k^iCki.tO^Oca.cO/^ki.k^q). (2.1) 

G is the colour factor for the process and the functions $j(ki,q) are the impact factors which 

contain the information about the external hadronic states. The function /(s,ki,k2,q) contains all 

the internal dynamics of the BFKL pomeron, is process independent and is well known in the leading 

logarithm approximation (LLA) which we use throughout [11, 12]. The corresponding differential 

cross-section is 

da _ \A(s,t)\ 2 



dt Wtts 2 



(2.2) 



For diffractive photon production we need to compute the 7*-7 impact factor. The relevant Feyn- 
man diagrams are shown in Fig. |l[ Dispersive techniques are used to calculate the imaginary part of 
the amplitude which is all that is required in the LLA. The study of the the photon impact factor is 



not new and dates back to [13|. In this paper, we compute the impact factor with an off-shell incoming 



photon and a general momentum transfer. 




Figure 1: The four diagrams that contribute to the photon-photon impact factor 



We detail the calculation of the contribution to the impact factor arising from the amplitude of 
Fig.|l](a). The same techniques apply to all diagrams. 
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The photon polarisation states are labelled by j and k and d(PS) refers to the two-body phase space 
of the cut quark lines shown in Fig.Q. We define two light-like and mutually orthogonal vectors, p\ 
and p2- All momenta are then Sudakov decomposed into components proportional to p%, and p2 and 
a third perpendicular vector: 
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(2.4) 

qi is the incoming photon momentum (Q 2 = —q 2 ), I the quark momentum, k\ + q and —k\ are 
the momenta of the gluon legs. The momentum transfer is q 2 = t < 0. Working in the eikonal 
approximation we are only interested in terms that are proportional to P\P\ in the numerator of (|2.3| ), 
Projecting out this piece, the amplitude becomes 
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(2.5) 



The two-body phase space integration over the cut lines is 
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(2.6) 



We have used the fact that p\ <C p in the LLA and have transfered from the Minkowski metric to the 
Euclidean one: 

- (^) 2 = k?- (2.7) 

The 6 functions are used to absorb the A and Ai integrals and lead to 

s 

ft = (Q 2 + q 2 )^ + ^. (2.8) 

s 



Within the LLA it is safe to neglect p\ and pt- 



Our notation follows that of [14|, and so we can readily extract the corresponding impact factor 
by removing the 
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where T 3 , -, is given by 
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It is not necessary to calculate the contributions to the impact factor from all four diagrams 
explicitly since they are related to each other by the transformations shown in Fig.||| 
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Figure 2: The relationship between the four amplitudes of Fig 



Thus, the full impact factor is given by 
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where 

jk tv(/ M/+ fa+ 4) tyjh 4i+ fa) My- 4i) 4)) 

w~ {i- qi y{i + k l + q y ' ( • } 

2.1 Photon polarisation 

The transverse polarisation vectors of the incoming photon lie in the transverse plane by definition. 
However, as the final state photon gets a kick into the transverse plane, its polarisation vectors have 
a component in the non-transverse directions. In the Lorentz gauge, the most general polarisation 
vectors for the incoming and outgoing photons are 



6 (/) ~ e ± - — P£> 

e% = e*/. (2.13) 

Placing these polarisation vectors into the amplitude and taking the trace [15| we get 



g } (ki,q) = ± 4a s a ^ e\ J dp j d 2 \ 2 
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(2.14) 



x = a,b,c,d labels the diagram in FigJIJ. The impact factor takes the + sign for diagrams (b) and (c) 
and the — sign for diagrams (a) and (d). We are focusing on diagram (c), in which case 

A( c ) = ki + pq. 

All boldface vectors are Euclidean and lie in the transverse plane, e.g. e 3 = — e^ . We have also 
introduced Q = Q 2 p(l — p). 

The longitudinal polarisation vector of the incoming photon is given by 

e = -J^Pi + P2 ■ ( 2 - 15 ) 



After substituting this expression into the amplitude and taking the trace: 

which is identically zero. Thus, in the LLA, longitudinally polarised photons do not scatter into real 
photons. 



2.2 Impact factor in position space 



The four-gluon amplitude, /(s,ki,k2,q) of (|2.1| ), can be written [12]: 

/(«,k 1 ,ka,q) = tAs / ^ 2 Pi ^P 2 d 2 p 3 d 2 p, e -i(k 1 -p 1 -(k 1+q ).p 2 -k 2 .p 3+ (k 2+q) .p 4 ) 
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r e^)^( Pl ,p 2 )^(p3,p 4 ). 



(2.17) 



(l/ 2 + l/4)2 

We have neglected the contributions with conformal spin n > since they give rise to sub-leading 
contributions at high energies. The eigenfunctions of the kernel are given explicitly by 

l+2iv 
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Xiy) = 2^(1) - tf (1/2 + iu) - ^(1/2 - iv) 
(^(x) is the digamma function) and 
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(2.18) 
(2.19) 

(2.20) 



The choice of scale sq is arbitrary in the LLA, in our subsequent calculations we choose either the 
largest scale (apart from s) in the problem or set a value for z. The number of colours, N c = 3. 
The scattering amplitude (|2.1| ) can now be written 



A(s,t) 



is G 



(27T) 



10 



dv 
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and we have introduced the functions 



11= /d 2 ^ ^i{\i X ^)Jd 2 Pl d 2 p 2 E»( Pl , P2 ) e -*X-(P 1 -p 2 )+iq-p 2 . 

For photon production we need to calculate ( 2.22J ) using &lZ- We work in the helicity basis: 

£ ± = --L ( i, ± i). 



There are four amplitudes to consider: (J, k) = (+, +), (+, —),(—,+),(—, — ) and 
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(2.21) 



(2.22) 



(2.23) 



(2.24) 



From now on, non-boldface momenta are understood to be complex numbers, e.g. I = l x + il y . To 
compute IY \ we first take the Fourier transform of the impact factor: 



${~f\p 1 ,p 2 )=e^Pi Jd 2 k, ^f\kx,q) e-^-(Pi-P 2 ). 
Summing over all four diagrams gives 
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Equivalently, and dropping for the time being the pre-factors and the integral over the momentum 
fraction p, we have 
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After some algebra, this becomes 
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(2.28) 



P + Q J \~ l J V + Q , 

where a and b are vector indices in the transverse space. Performing the angular parts of the 1 and 1' 
integrals and the ki and k 2 integrals yields 
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Using 



(2.29) 

(2.30) 
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Vr 2 J (|r 2 ||l|) = -|l| Ji(|r 2 ||l|)r 2 , 
e + ■ r 2 e~ ■ r 2 = 1, 
completing the radial 1 integral and re-arranging leads to the result: 

$(-,-)(Pi,P 2 ) = 16ir 2 aa s Ve 2 l l dp[p 2 + (l-pf\ /dWr 2 ^e^ ri e^-^ r2 

QKi(\r\Q) [<5 2 (n - Pl ) - 5\v 2 - Pl )\ [<5 2 (n - pa) - 5 2 (r 2 - p 2 )], (2.32) 
where r = ri — r 2 . The impact factors for the other polarisation states are 

$ (+,+)(Pl>P2) = ®(-,-)(Pl, P2) 

1 2 

$ (+) -)(Pi,P 2 ) = 327r 2 aa 8 X) e ?/ d P / d2r i d2r 2P( 1 -P)fi3 e * Pq ' Pl e* (1 - /,)q-Pa 
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$(-,+)(Pi,p 2 ) = 327r 2 aa s J24 f dp I d\ x d\ 2 p{\ - p) f^ e^' ri e^ 1 "'^ 

tfi(|r|Q) Q [J 2 (n - Pl ) - <5 2 (r 2 - Pl )\ [5 2 (n - p 2 ) - 5 2 (r 2 - p 2 )]. (2.33) 

3 Projection on Conformal Eigenstates 

We must now project the position space impact factors onto the conformal eigenfunctions of the BFKL 
kernel, i.e. we must perform the following convolution: 
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(3.34) 



This part of the calculation follows very closely that of ||]. For completeness, we reproduce the steps 
here. 



In two of the terms of ( |2.33D the ^-functions force p 1 = p 2 which results in a vanishing contribution 
due to the vanishing of E u . Hence, 
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(3.35) 



i.e. 
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q J ° 
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We next shift r 2 by — pr to eliminate the r dependent phase factor: 
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(3.37) 



We now switch to the complex notation, i.e from r = (77, r 2 ) to a = r\ + i r 2 and b = r\ — i r 2 , and 
use the freedom to choose q to be real: 
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(3.38) 



To proceed we factorise the integration in a, b and a 2 ,b 2 by writing the Bessel function as its Mellin 
transform, 
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and re-scaling a &; 6 by a 2 & b 2 respectively: 
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Re-scaling a 2 and b 2 by 2i/q gives 
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Completing the a, 6 and a 2 , o 2 integrals as in |B|: 
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and we have defined 
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The above integrals are representations of the hypergeometric function, i.e. ( p. 43 ) can be written 



2-F1 f-*2, *i + 1; zi + z 2 + 2; -^-J 2 Fi f-z 2 , -*i - 2z 2 - 1; -z\ - z 2 ; ^—j (3.45) 

and we have anticipated the symmetry between p and (1 — p) when performing the p integral. We can 
further simplify the hypergeometric functions to reach 
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Using these integrals, ( p. 41 ) reduces to 
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and we have used the identity 

r(z)T(l-z) = -^—. (3.48) 

The p-integration is performed by expanding the first hypergeometric function, performing the integral 
and then reducing the remaining series using $ (digamma) functions. The final result for the non-flip 
amplitude takes the following form: 
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and we have changed integration variables to z = 7 — 1 + i\. Of course the (— , — ) amplitude is 
equivalent. For the flip amplitudes we proceed in the same fashion. The final answer is 

T v (r fl 2x tv (r a 2n o- 5V 2{ 2 V~ 2w 1 tanh(vri/) 1 

J (+'-) (g '« } = J (-+) (0 '« } = "^^ ^ UJ r 2 (l/2 + ^) Tv (1 + , 2 ) 

/•7-l+ioo ( Q 2 \~ zl2 
/ ^py r(3/2 - »i/ - ar/2) r(3/2 + »i/ - */2) T(z/2) r(z/2 + l). (3.50) 

J7- 1— ioo I 9 / 

3.1 The limits Q 2 > g 2 and Q 2 < g 2 

In the limit Q 2 3> g 2 we are able to complete the z integration. We must close the contour in the 
right half plane. The leading contribution is extracted by considering the left most poles. 

For the non-flip amplitude these are located at z = 1 ± 2iv: 
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For the flip amplitude, the left most poles are located at z = 3 ± 2zz^ giving 
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We see that, in the limit Q 2 2> q 2 , the flip amplitude is suppressed by a factor of q 2 /Q 2 . 

In the limit Q 2 <C q 2 we must close the contour in the left hand plane. The leading piece is found 
by considering the pole at z = 0. This gives the Q 2 independent contribution found in ||]. The first 
correction to this is found by considering the double pole located at z = —2. We find 
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The analytic structure of the z plane is identical for the flip result and so 
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The cross section for 7*7* — > 77 scattering can be computed by substituting the results for I v given 
by ( 3.49] ) and ( p. 50 ) into ( 2.21 ). The colour factor for this process is given by 
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N 2 



2. 



(4.55) 



4.1 The small \t\ limit 

For this process we can study the |t| — > region. The cross section will be greatest in this limit and, 
interestingly, we can study the behaviour of the amplitude in \t\ — > limit. We focus on the limit 
Q 2 ^ Q 2 ^> \t\ where Q 2 and Q 2 , are the virtualities of the incoming photons. We shall find an analytic 
approximation for the 7*7* — > 77 amplitude using the approximation to I u given by equation ( |3.51| ) 
and evaluating the u integration using the saddle point method. We need only consider the non-flip 
amplitudes. When ( |3.51 ) is placed into (]2.21 ) we find two distinct terms. One term is independent 
of \t\ and gives the limiting value of the amplitude when |i| =0. The other term provides the leading 
|t| dependence. Considering the saddle point integration of the |t| dependent piece we find that the 
dominant term in the exponential is proportional to ln(Q 2 Q 2 /t 2 ) and so we may make the assumption 
that the saddle point lies at v ~ ±i/2 and approximate x( v ) by an expansion about v = ±i/2: 
x{ v ) ~ 1/(1/2 i iv)- The saddle point is in fact at 




4z 



ln(Q 2 Q 2 /t 2 ) 



± 1 -. 
2 



(4.56) 



For the \t\ independent piece we have a term proportional to h^Qf/Q 2 ,) i n the exponential. For 
Q 2 ~ Q\ the saddle point is at 



iHQl/Ql) 

28C(3)z 
10 



(4.57) 



and we can use xiy) ~ 4 log 2 — 14£(3)z^ 2 . Performing the integration over v yields 



A {+j+ . +i+) (s,t,QlQ 2 2 )\ Q2 ^ >>]ti wt 9* 2 a 2 a 2 a (5> 2 
/ vr 4 |t| In 5 / 4 (Q?Ql/|t | 2 ) 



Q1Q2 



exp 4zln2 



In 2 (gf/gj) 

56((3)z 



14C(3)z 9vr 7 /2 Q^a 



y 7/4 



exp 2Jzln(Q 2 Q 2 /|t 



2r>2/| + |2' 



(4.58) 



There is a cusp in the cross-section at —t = 0, i.e. the slope is infinite. In particular, defining an 

effective 6-parameter b e s by 

2 dA(s,t) 



Jeff 



A(s,t) dt 
and using ( [4.58 ), we find that, in the vicinity of —t = 0, 



1 4 1 



where 



A 4*ln2 / jj 



ln 5/4 ^OM) 



,7/4 



(4.59) 

(4.60) 
(4.61) 



4.2 The large \t\ limit 

We now consider the limit Q\-,Q\ <C \t\. The saddle point is at v = and, after integration over v, 
we find that 

/ \ 2 



A(+, + ; + ,+) (S, t, Ql,Q2, 



A (+t _. + _ ) (8,t,Qi,Qi 



\t\>Q\,Ql 



\t\»Q\,Ql 



1 s 

4 |t| VV 



22 '£e 2 i « 



4zln2 



/ 7T 



.3/2 



-7r— -a a„ 
4 t s 



E 



e 2 e 4 * ln2 



V14C(3)z 

/ vr \ 3 / 2 
U4C(3)z 



Ci(Qi)Ci(Q 2 ) 
C 2 (Qi)C 2 (Q 2 ) 



(4.62) 



where 



Cx(Qi] 

c 2 (Qi; 



11 + 



1Q| 

4 lil 



68 -307S -60 In 2 -15 In 



0? 



-lf(-f— «+- (I 



(4.63) 



The amplitudes A( +)+;+i _) and yl( + _. + + ) have a similar form, the final factors being replaced by 
Ci(Qi)C 2 (Q 2 ) and C 2 (Qi)Ci(Q 2 ) respectively. 



4.3 Numerical results 

In Fig. 3 we show the results for the 7*7* — > 77 differential cross-section obtained using equations 
( |3.49| ) and fl3.50|) and compare with the results obtained using the approximations to IY, ,-, and IY, % 
given in (|3.51| ), ( [3.52 ), ( 3. 53| ) and (3.54). All v integrals are computed numerically. 

At LEP2, z values around 1 are typical and values of around 2 are more typical of those attainable 
at a future linear collider. 
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5 YP -► 7^ 



We may use the results ( p.49|) and (|3.50|) to compute 7* —>■ 7 scattering off a hadronic state. We 
consider 7*P — ► 7X scattering at the energies typical of those at the HERA collider. 

The impact factor for the coupling to the quarks and gluons in the proton is well known j^, [uj] 
and yields 

, x = T(l/2-iu) /2\ 1 ~ 2iu 

5 = -w "-WftTiJj (J (5 ' 64) 

for the coupling to quarks. The colour factor for coupling to quarks is 

N 2 -l 

relative to which the coupling to gluons is enhanced by a factor Ca/Cf = 9/4. Putting the colour 
factor, along with (|3l9|),(p0|),([O4D, into ( pip yields 



A 



(+,+) 



[s,t,Q 2 



i aa 2 y^ e 2 



7T s 

9 6 W\ 



du 



1 tanh(7rz/) 



(1/4 + 1/2)2 1 + zy 2 ^ 



,«xM 



7-l+ioo dz /Q2- 
7— 1— xoo ^7T2 \ |£| 



-z/2 



r(l/2 -iv- z/2) r(l/2 + ti/ - z/2) 
r(l/2 + ii/) r(l/2-ii/) 



r(z/2) r(z/2 + l)[z 2 + 11 + 12i/ 2 ] 

(5.65) 



and 



A {+ _) (s, i, Q 2 ) = i aa 2 s Y,e 2 q Y W\ j dl/ 



1 tanh(7rz/) 



^xM 



7 -l+*oo dz /q2~ 
7— 1— zoo ^7TZ \ \t\ 



q 

z/2 



(1/4 + u 2 ) 2 l + i/ 



tt;/ 



r(3/2 - fi/ - z/2) r(3/2 + ti/ - z/2) 
r(l/2 + ii/) r(l/2-ii/) 



T(z/2) T(z/2 + l). (5.66) 



5.1 The low Q 2 limit 

In the limit —t ^$> Q 2 the saddle point integration about v = yields 



A (+i+) (s,i,Q 2 



87 



/ vr \ 3 / 2 



"^Ve^e 42 ")— t] Ci(g) 
? 2 «l<l 3 „ q \t\ V14zC(3)/ WJ 



(5.67) 



and 



^ (+ ,_)(s,t,Q 2 



8i , x-^ 9 s 4z In 2 / 7r 



2 «l*l 3 



ot 2 x - ^ 2 ° 

— aa„ > e„ 7-7 e' 



14zC(3) 



3/2 



C 2 (Q). 



(5.68) 
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5.2 The high Q 2 limit 

The relevant saddle point is now at v ~ ±i/2 and so 



8t 2 _ 2 s ln 3 / 4 (Q 2 /|i|) 



A (+ , +) ( S ,t,Q*) Q2>|t| = 3^= ao^ej qs ^71^ ex P [^ *HQ 2 l\t\)) (5.69) 

and 

A+,-)M,Q 2 )| Q2>>|t| = ^ ««?E^ ^ ^gyi^ exp (2^1n(gVI*l)) (5-70) 

The dominance of the non-flip amplitude over the flip amplitude is clear: 

A (+,+) _ 3Q 2 



A(+,-) 2|t| 



(5-71) 



5.3 Numerical Results 



In Fig. 4 we show the results for proton-photon scattering in by performing the v and z integrals 
exactly. We compare with the results obtained using the analytic approximations to the amplitude 
derived from (|3.52j) ,( J3~l)l" ), ( 3.541 ), (|3.53| ) again with the v integral performed exactly. 



The proton-photon cross section is determined by convoluting the photon-quark cross section with 
the parton density functions [17, 18|: 



dajjP -> jX) 
dt 



£ dx fc±g(x, t) + J2(q( X , t) + q(x, t))) ^ (7 y 7g) . (5.72) 



The lower limit of the x integral is determined by the on-shell condition for the struck quark in the 
proton, i.e. 

' -, (5.73) 






where Mx is the maximum invariant mass of the proton dissociation products. We choose Mx = 10 
GeV, which is typical of the cuts imposed by the HERA experiments. We note that in the limit of 
|t| S> Q 2 the cross section becomes insensitive to values of Q 2 , as anticipated in the analytic expressions 
( |5.67 ) and ( |5.68| ). The approximate expressions are shown to be good over a large range in Q 2 /\t\. 



We note that these cross-sections are large enough to make observation at HERA feasible. 
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Figure 3. 7*7* — ► 77 cross section comparing the exact BFKL result with the leading 
approximation. The solid line gives the exact result and the dotted line the result obtained by only 

keeping the leading terms as described in the text. 
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7-P cross section BFKL LO Exact. G 2 =0,' ,10,100 GeV 2 ,a, = 0.2 
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Figure 4. 7*P — ► jX cross section comparing the exact BFKL result with the leading 
approximation. The solid line gives the exact result and the dotted line the result obtained by only- 
keeping the leading terms as described in the text. 
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